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Abstract: The lightest neutralino of -R-parity conserving supersymmetric models serves 
as a compelling candidate to account for the presence of cold dark matter in the universe. 
In the minimal supergravity (mSUGRA) model, a relic density can be found in accord with 
recent WMAP data for large values of the parameter tan /3, where neutralino annihilation in 
the early universe occurs via the broad s-channel resonance of the pseudoscalar Higgs boson 
A. We map out rates for indirect detection of neutralinos via 1. detection of neutrinos 
arising from neutralino annihilation in the core of the earth or sun and 2. detection of 
gamma rays, antiprotons and positrons arising from neutralino annihilation in the galactic 
halo. If indeed ^4-resonance annihilation is the main sink for neutralinos in the early 
universe, then signals may occur in the gamma ray, antiproton and positron channels, 
while a signal in the neutrino channel would likely be absent. This is in contrast to the 
hyperbolic branch/focus point (HB/FP) region where all indirect detection signals are 
likely to occur, and also in contrast to the stau co-annihilation region, where none of the 
indirect signals are likely to occur. 
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In the past decades, a growing body of astrophysical evidence has made an irrefragable 
case for the existence of cold dark matter (CDM) in the universe The most recent results 
come from the Wilkinson Microwave Anisotropy Probe (WMAP)|Q]. Their results confirm 
the standard model of cosmology and fit its parameters to high precision. The properties 
of a flat universe in the AC DM model are characterized by the density of baryons (fib), 
matter density (fi m ), vacuum energy and the expansion rate (h) which are measured 
to be: 

Q b = 0.044 ± 0.004 (1) 
n m = 0.27 ± 0.04 (2) 
n A = 0.73 ± 0.04 (3) 

h = o.7ij5:8*. W 

From the WMAP results, a value for the cold dark matter density of the universe can be 
derived: 

ncDMh 2 = o.ii26t8:889o(^o:oi8i) a * 68(95)% cl. (5) 

A particularly attractive candidate for CDM is the lightest neutralino in /2-parity 
conserving supersymmetric models Q. In the paradigm minimal supergravity (mSUGRA) 
model P], it is assumed that at the scale Q = Mqut, there is a common scalar mass 
mo, a common gaugino mass mi/2, and a common trilinear term Aq. The soft SUSY 
breaking terms can be calculated at scale Q = M wea k via renormalzation group evolution. 
Electroweak symmetry breaking occurs radiatively (REWSB) due to the large top quark 
mass, so that the bilinear soft breaking term B can be traded for the weak scale ratio of 
Higgs vevs tan (3, and the magnitude (but not the sign) of the superpotential [i term can 
be specified. Thus, the mSUGRA model is characterized by four parameters plus a sign 
choice: 

mo, ^1/2) A); tan/3, and sign(n). (6) 

Once these model parameters are specified, then all sparticle masses and mixings are de- 
termined, and scattering cross sections may be reliably calculated. 

In the early universe at very high temperatures, the lightest neutralino Z\ will be in 
thermal equilibrium, so that its number density is well determined. As the universe expands 
and cools, there will be insufficient thermal energy to produce neutralinos, although they 
can still annihilate with one another. The neutralino relic density can be determined by 
solving the Boltzmann equation for neutralinos in a Friedmann- Robertson- Walker universe. 
In spite of the claims that the lightest neutralino is a good dark matter candidate, it turns 
out that in most of the parameter space of the mSUGRA model, a value of f^^i 2 well 
beyond the WMAP bound is generated. Only certain regions of the mSUGRA model 
parameter space give rise to a relatively low value of £1% h 2 in accord with astrophysical 
measurements and theory. These regions consist of : 

1 Additional less prominent parameter space regions are also possible, such as the light higgs h resonance 
region (bordering the LEP2 bounds at low m 1 f 2 ) and the top squark co-annihilation region (for very 
particular choices of the Aq parameter). 



- 1 - 



1. The bulk annihilation region at low values of mo and mi/2> where neutralino pair 
annihilation occurs at a large rate via ^-channel slepton exchange. 

2. The stau co-annihilation region at low mo where m^ ~ m^ so that Z±s may co- 
annihilate with fis in the early universe ||. 

3. The hyperbolic branch/focus point (HB/FP) regionQ at large mo near the boundary 
of the REWSB excluded region where \/j,\ becomes small, and the neutralinos have 
a significant higgsino component, which facilitates annihilations to WW and ZZ 
pairs p. 

4. The ^-annihilation funnel, which occurs at very large tan/3 ~ 45 — 60|7|. In this 
case, the value of tua ~ 2m ^ . An exact equality of the mass relation isn't necessary, 
since the A width can be quite large (r^ ~ 10 — 50 GeV); then 2m^ can be several 
widths away from resonance, and still achieve a large Z\Z\ —> A —> ff annihilation 
cross section. The heavy scalar Higgs H also contributes to the annihilation cross 
section. 



Several years ago, the bulk annihilation region of parameter space was favored. This 
situation has changed in that the bulk annihilation region generally predicts a light SUSY 
Higgs boson h with mass below LEP2 bounds, along with large- usually anomalous- pre- 
dictions of the rate for BF(b —* sj) decays and muon anomalous magnetic moment 
o-fj, = id — 2)^/2 ]s|, [|. An increase of either of the parameters mo or m\ii leads generally 
to heavier sparticle masses and m^ values, so that predictions for loop induced processes 
become more SM-like. 

Given a knowledge of which regions of model parameter space give rise to neutralino 
relic densities in accord with measurements, it is useful to examine the implications for 
detection of supersymmetric matter. The direct sparticle search limits for the Fermilab 



Tevatron collider |L0|], the CERN LHC[[Ll| and a y/5 = 0.5 - 1 TeV linear collider^] have 
all been examined. In addition, there exist both direct and indirect dark matter search 
experiments that are ongoing and proposed. Direct dark matter detection has been recently 



examined by many authors] 13 1, and observable signal rates are generally found in either 
the bulk annihilation region, or in the HB/FP region, while direct detection of DM seems 
unlikely in the j4-funnel or in the stau co-annihilation region. 

Indirect detection of neutralino dark matter may occur via 

1. observation of high energy neutrinos originating from ZiZi annihilations in the core 
of the sun or earth [15], 



2. observation of 7-rays originating from neutralino annihilation in the galactic core or 
halo pi] and 



3. observation of positrons [|l7j or anti-protons [18] originating from neutralino annihila- 
tion in the galactic halo. 
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The latter signals would typically be non-directional due to the influence of galactic mag- 
netic fields, unless the neutralino annihilations occur relatively close to earth in regions of 
clumpy dark matter. 

The indirect signals for SUSY dark matter have been investigated in a large number 
of papers, and computer codes which yield the various signal rates are availablepl, M. 
Recent works find that the various indirect signals occur at large rates in the now disfavored 
bulk annihilation region, and also in the HB/FP region[21]. Naively, this is not surprising 
since the same regions of parameter space that include large neutralino annihilation cross 
sections in the early universe should give large annihilation cross sections as sources of 
indirect signals for SUSY dark matter. 

In this paper, we pay special attention to indirect signals for SUSY dark matter in the 
A-annihilation funnel. We generate sparticle mass spectra using Isajet v7.69[22], which 
includes full one-loop radiative corrections to all sparticle masses and Yukawa couplings, 
and minimizes the scalar potential using the renormalization group improved 1-loop effec- 
tive potential including all tadpole contributions, evaluated at an optimized scale choice 
which accounts for leading two loop terms. Good agreement between m^ values is found 
in comparison with the FeynHiggs program, and there is good agreement as well in the 
rriA calculation between Isajet and SoftSUSY, Spheno and Suspect codes, as detailed in 
Ref. 22]. To evaluate the indirect signals expected from the mSUGRA model, we adopt 
the DarkSUSY package[20] interfaced to Isasugra 2 . For our calculation of the neutralino 
relic density, we use the Isared program 1 24] interfaced with Isajet. Isared calculates all 
relevant neutralino pair annihilation and co-annihilation processes with relativistic ther- 
mal aver aging f25j| . An important element of the calculation is that Isared calculates the 
neutralino relic density using the Isajet 2-loop t, b and r Yukawa couplings evaluated at 
the scale Q = itia- The Yukawa coupling calculation begins with the DR fermion masses 
at scale Q = Mz, and evolves via SM renormalization group equations (RGEs) to the scale 
Qsusy = \J m i L m t R i where complete MSSM 1-loop threshold corrections are implemented. 
Evolution at higher mass scales is implemented via 2-loop MSSM RGEs. The final RGE 
solution is gained after iterative running of couplings and soft terms between Mz and 
Mqut and back until a convergent solution is achieved. 

We first show in Fig. [l| the location of the ^4-pole in the mSUGRA model by plotting 
\niA — 2m^ \/Ta versus mo for mi/ 2 = 500, 750, 1000 and 1250 GeV, for Aq = and a) 
tan (5 = 45 and /U < and b) tan (3 = 54 and /i > 0. The point in tuq where \rriA — ^ m z 1 \/^A 
drops to zero shows the location of the A-annihilation funnel. One may also note from Fig. 
|l] that there is a band of mo values wherein 2m^ is within several widths of the A-pole. 
As tan (3 increases, the b Yukawa coupling also increases, leading to large values of the A 
width r^. For very large values of tan/3 ~ 54, Ta can exceed 50 GeV. In frame a), it can 
be seen that the ^-annihilation funnel occurs for all values of m^ 2 shown. However, in 
frame 6), it can be seen that the ^4-annihilation resonance is only found for large values 
of mi/2 ~ 1 TeV. Even so, for lower m^j values, the effect of the ^-annihilation funnel is 
felt at the low mo range since one may be only one to several widths away from resonance. 

2 Isasugra is a subprogram of the Isajet package that calculates sparticle mass spectra and branching 
fractions for a variety of supersymmetric models 
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The explicit location of the ^4-annihilation funnel in the rriQ vs. plane as calculated 

by Isasugra is shown in Ref. 0, and will not be repeated here. 
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Figure 1: A plot of \m,A — SmgJ/r 'a versus mo for various mi/2 values and Aq — and a) 
tan (3 = 45 with /i < and b) tan/3 = 54 and /i > 0. 

In Fig. [2|, we show in the top frames the neutralino relic density O ^ h 2 versus mSUGRA 
parameter mo for various m^ 2 values. We also take Aq = and a) tan/3 = 45 and /i < 
and b) tan/5 = 54 and fi > 0. In frame a), we see that most of the range of mo yields a 
value of f2 g h 2 far beyond the WMAP measured range (shown by horizontal dotted lines) . 
However, as mo becomes very large, one enters the HB/FP region, and |^| becomes small. 
The Z\ gains a substantial higgsino component which facilitates Z\Z\ annihilation into final 
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states such as WW, ZZ and Zh. Thus, ^j^^ 2 drops into and below the WMAP measured 
range for CDM. As too decreases to small values, there is also a decrease in neutralino relic 
density. This time the decrease is due to the fact that to ^4 decreases towards the value 
2m^, so that Z\Z\ — > A, H —* bb is enhanced. For low m\ji values, Q^^ 2 actually drops 
to very low values — below 0.025 — which wouldn't even be enough to explain galactic 
rotation curves. However, the A-annihilation funnel is quite broad because the b and r 
Yukawa couplings grow with tan /3, which increases the A and H widths, and also the Z\Z\ 
annihilation rate through s-channel A and H exchange. Optimal values of ^^_h 2 in accord 
with the WMAP CDM result are then achieved close to resonance, but not exactly on it. 
For even lower values of mo, the value of ^^ ± h 2 increases as one moves away from the A 
resonance, until there is a final down-turn in Q^Jt? due to stau co-annihilation. In frame 
b) at the top, we show the relic density for \x > and tan/3 = 54. Again, much of the range 
in too is excluded since Q^h 2 is beyond the WMAP limit. For large too, the relic density 
again drops as a HB/FP region is approached. For low values of too, the relic density 
again decreases due to the growing importance of neutralino annihilation via s-channel A 
and H exchange. However, only for mx/2 ~ 1 TeV do we actually meet the A resonance. 
Nonetheless, for lower mx/2 values, the effect of the A and H pole is felt, and decreases the 
relic density to sub- WMAP values until the stau co- annihilation region is hit at the very 
lowest values of too- 

In frames c) and d), we show the neutralino-proton spin independent (scalar) scattering 
cross section asi for tan/3 = 45, \x < and tan/3 = 54, /i > 0, respectively. For low values 
of too, the scattering cross section is enhanced because squark masses become relatively 
light, and u-channel squark exchange graphs yield large scattering amplitudes. For high 
too values, the neutralino has a significant higgsino component, which enhances the t- 
channel Higgs exchange diagrams. For the negative \x case in frame c), there is destructive 
interference between Higgs and squark exchange diagrams, and the cross section drops to 
zero for a particular too value. In frame d), there is no destructive interference, so the 
spin-independent cross section just drops to a minimal but finite value. We note here that 
Stage 3 dark matter detectors p3[ such as Cryoarray, Zeplin-4 and Genius aim towards a 



sensitivity of roughly 10~ 9 pb, depending somewhat on the value of m^ . We also note 
that the neutralino accretion rate for the earth depends strongly on the spin-independent 
neutralino-nucleon scattering, which is enhanced for the heavy nuclei of which the earth is 
composed. 

In frames e) and f), we show the neutralino-proton spin-dependent (axial- vector) scat- 
tering cross section asD versus too for the same mSUGRA parameters as the previous 
figures. The cross section is again enhanced in the HB/FP region, but drops to a mini- 
mum as too decreases. We note here that the neutralino accretion rate of the sun depends 
strongly on the spin-dependent neutralino-nucleon scattering cross section. 

In Fig. ^, we show the flux of muons coming from neutralino annihilation in the 
core of the earth (frames a) and b)), and from neutralino annihilation in the core of the 
sun (frames c) and d)), for the same mSUGRA parameter values as in Fig. Second 
generation neutrino telescopes such as Icecube(2^] and Antares|27] hope to probe muon flux 
values of 10-100 muons/km 2 /yr|L4|]. Comparing this number to the results from frames a) 



- 5 - 



10' 

10' 

10° 
t\l , 

.c 10 

a 

10" 
10^ 
10" 4 

10" 6 

io- 8 



Q. 

J"10" 1 



10"* 

io- 16 

10" : 
10" 



Q. 

I10- 5 



10" 

io- ; 



tanp =45 

\i<0 





1 i 1 i 1 i 




\ Ef 


V \ \ \ i 


i 


1,1,1," 



2000 



2000 



4000 



1 1 


1 1 1 
/C^Z^^' m„ 2 =300 GeV - 




— m 1I2 =400 GeV 




m, a =500 GeV 




m 1:2 =600 GeV 




m, e =750 GeV 




m 1I2 =1000 GeV - 


c) 'I 

1 


— m 1(2 =1250 GeV 





4000 




2000 4000 

m (GeV) 




6000 



6000 



2000 4000 

m (GeV) 



6000 



Figure 2: We show the neutralino relic density versus mo for various values of fn\ji and a) 
tan/3 = 45, fi < and 6) tan/3 = 53, \i > 0. We take Ao = 0. In frames c) and d), we show 
the spin-independent neutralino-proton scattering cross section, while in frames e) and f), we show 
the spin-independent neutralino-proton scattering cross section versus too for the same mSUGRA 
parameters as used in frames a) and b). 



and b) shows that a signal from the mSUGRA model is unlikely to come from neutralino 
annihilation at the center of the earth. The drop in muon flux at moderate mo values follows 
along the curves shown previously for the neutralino-proton spin-dependent scattering cross 
section. There is an enhancement in muon flux at low and high mo values, but probably 
not enough to create a detectable signal. We note here that the muon flux from the earth 
shows some enhancement in rate in the ^-annihilation funnel, but not enough to push the 
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expected flux levels into the observable regime. The solid curves are plotted assuming a 
local relic density of neutralinos given by p$ = 0.3 GeV/cm 3 . If the relic density falls to 
very low values, then the local density may have to be rescaled in accord with the global 
relic density. The dashed curves show the variation in the rates if the local relic density is 

rescaled according to p = po Q Q 1 25 for Vl^ h values below 0.025. 
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Figure 3: Muon flux from the earth (frames a) and b)) and sun (frames c) and gQ) versus mo for 
various values of m 1 / 2 and tan/3 = 45, p < (left-hand frames) and tan/3 = 53, p > (right-hand 
frames). We take Aq = 0. Dashed lines include rescaling of the local relic density for values of 
tt% h 2 < 0.025. 

Alternatively, if we examine frames c) and d), we see that detectable levels of muon 
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flux may indeed arise from neutralino annihilations at the core of the sun. These large rates 
occur in the HB/FP region, as noted by other authors [21]. It is intriguing that the large 
muon flux occurs in one of the main regions where the relic density is in accord with WMAP 
analyses. However, in the ^-annihilation funnel, there is no evidence of enhancement of 
the muon flux. This is because the rate for neutralino annihilation in the sun or earth is 
given by 

T A = ^Ctanh 2 (VCAt Q ), (7) 

where C is the capture rate, A is the total annihilation rate times relative velocity per 
volume, and i is the present age of the solar system. For the sun, the age of the solar 
system exceeds the equilibration time, so Ta ~ and the muon flux tends to follow the 
neutralino-nucleon scattering rate rather than the neutralino pair annihilation cross section. 
The earth has typically a much longer equilibration time, so that Ta ~ ^C 2 At 2 , and is 
hence more sensitive to the neutralino annihilation cross section times relative velocity. 

In Fig. |], we show rates for gamma rays (frames a) and b)), positrons (frames c) 
and e?))and anti-protons (frames e) and J)) originating from neutralino annihilation in the 
galactic core and halo. The plots are shown versus mo for the same mSUGRA parameters 
as in Figs. [2] and pL In frames a) and 6), the flux for continuum gamma rays with 
energy 2£ 7 > 1 GeV is shown in units of photons /cm 2 /sec, assuming a detector with 0.001 
sr solid angle coverage, pointed at the galactic center. For halo model dependence in 
the distribution of dark matter, we adopt default DarkSUSY values. Experiments such 
as GLAST[28| expect to probe flux rates as low as 10 -10 photons/cm 2 /sec. Thus, in 
Fig. ||, we see that observable rates are expected to occur in the HB/FP region for both 
tan j3 = 45 and tan (3 = 54 cases. In addition, observable rates are expected if SUSY model 
parameters lie in the ^-annihilation funnel. At low mo values, the gamma ray flux rises 
and follows the ^-annihilation resonance. Detectable rates may occur in the j4-funnel as 
long as my2 ~ 1 TeV. We also show again as dashed curves the rates if the relic density is 
rescaled when ^^_h 2 falls below 0.025. In this case, rates for gamma ray detection may fall 
below observable levels, but only because neutralinos would not be the major constituent 
of CDM. If ^zji 2 lies within the WMAP band, then typically 2m^ will lie somewhat 
off-resonance, and observable rates for gamma ray detection can be found. In frame b) for 
fi > 0, detectable rates again occur at low mo for m l / 2 values below 500 GeV. 

In Fig. |4] frames c) and d) we show the signal-to-background (S/B) rates for detection 
of positrons arising from neutralino annihilations in the galactic halo. To calculate the 
S/B rates, we adopt fit C from Ref. Q for the E 2 d<fr e + /dttdE background rate. We 
compute the signal using the DarkSUSY positron flux evaluated at an "optimized" energy 
oiE = /2, as suggested in Ref. @. A S/B ~ 0.01 rate may be detectable pi], by 
experiments such as PAMELA |29| and AMS-02pO[. We see from Fig. | that observable 
rates may again occur in the HB/FP region, and also in the yl-annihilation funnel. 

Finally, in frames e) and f), we show the differential flux of antiprotons from the 
galactic halo, d&p/dEpdQ, for Ep = 1.76 GeV. In this case, background rates are more 
uncertain, so we show only the differential flux of antiprotons/GeV/cm 2 /sec/sr. Again, 
the largest rates occur in the HB/FP region, and also in the ^-annihilation funnel. 
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Figure 4: Flux of continuum gamma rays from a O.OOlsr cone centered on the galactic center, 
with Ey > 1 GeV (frames a) and 6)). We also show S/B for positrons (frames c) and dj) and the 
differential flux of antiprotons with Ep = 1.76 GeV (frames e) and f)). All plots are versus mo for 
various values of ra\ii and tan/3 = 45, \i < (left-hand frames) and tan/3 = 53, /i > (right-hand 
frames). We take Aq — 0. Dashed lines include rescaling of the neutralino density for values of 
Q,y h 2 < 0.025. 



We can summarize our results according to mSUGRA parameter space regions which 
give rise to a reasonable relic density of CDM. In the HB/FP region, both direct and 
indirect detection of SUSY dark matter is possible. Indirect detection signals may be 
possible for neutrino telescopes by observing muons arising from high energy neutrinos 
which are produced from the decays of final state originating from neutralino annihilation in 
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the core of the sun. Detectable rates also may occur for cosmic gamma rays, positrons and 
possibly anti-protons. If instead the A-annihilation funnel is the main sink for neutralinos 
in the early universe, then it is unlikely any signal would be seen from high energy neutrinos 
arising from neutralino annihilation in the core of the sun or the earth. However, observable 
signals in the gamma ray, positron and possibly antiproton channels may occur. Finally, if 
the sink for early universe neutralinos is due to stau co-annihilation, then there may be no 
direct or indirect signals for SUSY DM. The exception occurs at very large tan (3 values, 
where the ^4-annihilation funnel and stau co-annihilation region begin to overlap, in which 
case gammas, positrons and antiprotons could be visible, while neutrino-induced muons 
would not be. 
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